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(I)  Pm^ODUCTION 


The  research  program  undertaken  for  this  ■Report  is  divided  into  three 
parts.  In  the  first  part  we  analyze  the  simultaneous  observations  of 
different  mesospheric  optical  emissions  within  the  ongoing  MAPSTAR  program 
in  which  we  have  also  had  seme  participation  in  the  field  work.  We  then  in 
the  second  part  decompose  these  observations  into  their  monochromatic 
Fourier  components  and  match  them  with  COSPAR-based  gravity  wave  models  to 
obtain  various  wind  profiles  from  which  other  atmospheric  parameters  (e.g. 
the  kinetic  energy  density)  can  be  derived.  The  results  can  be  compared 
with  the  more  direct  simultaneous  observations  of  wind  profiles  from  MST 
radar,  which  will  form  part  of  the  future  MAPSTAR  program.  Such  comparisons 
can  connect  the  ground-based  optical  emission  data  with  the  simultaneously 
observed  radar  data  and  provide  information  on  gravity  wave  models  as  well 
as  on  the  chemistry  and  dynamics  of  wave- induced  airglcw  fluctuations.  In 
this  report,  we  shall  present  the  steps  which  can  allcw  us  to  extract  the 
velocity  fields  from  the  observed  airglcw  using  the  purely  dynamical  model 
of  Hatfield  et  al  (1981) .  A  proper  chemical -dynamical  model  is  under 
development  and  the  methods  presented  here  can  be  immediately  adapted  to 
the  more  fully  developed  model. 

In  part  three  we  begin  an  investigation  into  the  behavior  of  small 
scale  spatial-temporal  airglcw  fluctuations  induced  by  gravity  waves.  The 
investigation  is  of  intrinsic  physical  interest  in  that  it  examines  a 
simultaneous  ducting  and  filtering  mechanism  for  small-period  gravity  waves 
which  are  restricted  to  lewer  height  levels,  leaving  the  longer  period 


waves  to  penetrate  to  higher  altitudes.  In  terms  of  airglow,  such  a  result 
can  imply  an  overall  greater  probability  of  finding  small  scale  spatial- 
temporal  mesospheric  airglcw  fluctuations  occurring  simultaneously  with 
larger  scale  ionospheric  airglcw  fluctuations.  This  applies  only  to  lower 
altitude  sources  and  some  of  the  medium-scale  T.I.D. *s.  The  problem  should 
also  be  of  interest  to  the  Air  Force,  since  it  deals,  in  particular,  with 
snail  scale  infrared  fluctuations. 

For  the  first  part,  we  shall  deal  specifically  with  the  observations 
and  analyses  of  the  7300A  OH  bands,  the  5577A  01  green  line,  as  well  as 
some  radiometric  data  with  a  1.53  ^m  window.  The  observations  have  been 
made  at  Colorado  (Bootlake) ,  Wyoming  (Jelm)  and  New  Mexico  (Capilla  Peak) 
by  Peterson,  Pendleton  and  Lin,  as  part  of  the  MARSTAR  canpaign  of  1984  and 
most  of  the  observations  were  taken  during  a  period  from  May  20  to  July  3, 
1984.  The  reason  is  that  for  only  these  two  particular  days  we  have  both 
the  radiometric  and  the  photometric  data  as  well  as  even  some  interesting 
interferometric  data  which  sure  unfortunately  not  so  easily  usable  because 
of  the  relatively  short  and  interrupted  data  base.  There  is  also  the 
additional  reason  that  day  176  represents  a  typically  active  day  while  day 
182  represents  a  very  active  day.  Actually,  both  from  the  observed 
fluctuations  and  from  the  computation  of  percentage  fluctuations  day  180 
(June  27,  1984)  is  unquestionably  the  most  disturbed  day  of  the  1984 
canpaign,  but  since  we  do  not  have  the  radiometric  data  for  day  180  for  a 
more  complete  analyses,  we  shall  work  with  day  182  as  the  next  best 
exanple. 
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Ttie  analyses  will  be  divided  into  three  categories: 

(1)  Fhase  Correlation 

(2)  Frequency  Correlation 

(3)  Comparison  of  percentage  r.m.s.  fluctuations 

For  (1)  we  shall  superimpose  photometric  and  radiometric  observations 
made  during  the  same  time  period  and  at  the  same  location  to  establish  any 
phase  correlation  between  different  spectral  lines.  Cross  covariances  will 
also  be  computed  to  lend  more  quantitative  support  to  the  simple  visual 
comparisons.  As  we  shall  see,  the  phase  relationships  can  be  very 
complicated  depending  on  factors  ranging  from  the  nature  and  location  of 
the  sources,  the  geometry  of  the  observation  point  relative  to  the  source, 
to  the  shape  of  the  profiles  of  the  ambient  atmospheric  constituents  which 
produce  these  optical  emissions  (Theme  (1968),  Porter  &  Tuan  (1974),  Chiu  & 
Ching  (1978)  &  Hatfield  et  al  (1981)). 


For  (2)  we  shall  use  four  different  kinds  of  power  spectrum  analyses; 
the  maximum  entropy  (M.E. ) ;  the  Blackman  and  Tukey  (B.T. ) ;  the  modified 
periodograms  and  the  M-Step  Difference  Filter  can  M.E.  we  shall  base  our 
analysis  primarily  on  the  M.E.  and  use  the  others  merely  as  a  check  since 
the  M.E.  has  the  best  resolution  for  the  size  of  our  data  base  with  the 
Nyquist  frequency  corresponding  to  a  period  of  one  minute  which  is  well 
below  the  smallest  period  we  seek  (about  5  min.) .  The  results  of  the 
analyses  will  then  be  listed  (after  subtracting  the  Window  affects)  with 


the  most  prominent  periodicity  underlined.  We  shall  also  list  the  common 
periodicity  between  different  spectral  lines  observed  at  the  same  location 
and  time.  Unlike  the  comparison  of  phase  which  depends  so  much  on  the 
shape  of  the  ambient  undisturbed  background  intensity  distribution,  the 
frequencies  are  independent  of  the  shape  of  the  ambient  background  and 
common  frequencies  may  be  expected  to  occur  even  if  there  does  not  appear 
to  be  a  definite  phase  relationship. 

For  (3)  we  shall  list  the  percentage  r.m.s.  fluctuations  in  the  air 
glcw  averaged  over  60,  30,  15,  10  and  5  minutes.  Since  there  is  no  good 
unique  method  to  determine  the  general  trending  which  may  occur  for  reasons 
other  than  purely  atmospheric  dynamics,  we  have  decided  to  average  over  the 
above  sets  of  different  time  intervals.  The  average  over  the  shorter 
periods  will  automatically  include  trending  without  the  necessity  to 
determine  the  cause  for  such  trending.  The  averaging  over  different  time 
intervals  will  also  allcw  us  to  evaluate  the  relative  importance  of  small 
period  oscillations  as  compared  to  the  longer  period  oscillations.  For 
instance,  the  long  period  fluctuations  of  one  day  may  be  more  significant 
than  another,  but  the  short  period  fluctuations  may  be  precisely  the 
reverse.  The  comparison  of  r.m.s.  percentage  fluctuations  will  also 
provide  us  with  a  more  quantitative  means  of  picking  out  the  most  active 
day  as  well  as  the  least  active  day.  Lastly,  just  as  the  phase 
relationship  between  different  spectral  lines  depend  on  the  shape  of  the 
ambient  profile,  the  magnitude  of  the  response  of  the  various  "minor" 
atmospheric  constituents  (which  produce  the  atmospheric  emission)  to  an 
incoming  gravity  wave  also  depends  on  the  shape  of  the  ambient  density 

If 


A 


profile.  In  general,  the  response  is  much  greater  if  the  vertical  density 
gradients  is  positive,  i.e.  opposite  to  the  vertical  density  gradients  of 
the  major  atmospheric  constituents.  Thus,  fluctuations  in  the  optical 
emission  may  well  be  greater  owing  to  the  peculiarities  of  its  ambient 
emission  profile.  If  all  or  most  of  the  emission  profile  has  a  layered 
structure  with  a  very  narrow  half-width,  then  the  very  large  vertical 
gradients  would  produce  a  rather  large  fluctuation.  This  means  that  the 
percentage  r.m.s.  fluctuation  would  also  tend  to  be  more  significant. 


The  second  part  deals  with  the  extraction  of  velocity  field  profiles 
from  the  observed  data.  We  have  decomposed  the  observed  airglcw  data  for 
June  23  and  June  29,  1984  into  their  Fourier  components  from  FFT  power 
spectrum  analyses.  For  the  present  case  we  cannot  use  the  Max.  entropy 
method  for  the  simple  reason  that  it  does  not  lend  itself  readily  for 
Fourier  decomposition.  We  will  instead  use  the  feist  Fourier  with  a  Haming 
window.  Since  the  undisturbed  background  emission  is  a  very  slowly  varying 
function  of  time,  the  background  will  be  approximated  by  its  time-averaged 
value.  In  this  way  we  can  obtain  the  percentage  fluctuation  for  each 
Fourier  component  corresponding  to  a  particular  frequency.  We  can  then  use 
the  theoretical  computations  of  the  fluctuations  of  airglcw  (Hatfield  et  al 
(1981))  (Modified  to  include  the  fluctuation  of  the  chemical  rate 
coefficient  from  temperature  waves)  produced  by  CDSPAR  gravity-wave  models 
for  the  same  fixed  frequency  and  adjust  the  amplitudes  of  the  vertical  and 
horizontal  velocity  wave  components  until  the  theoretically  computed 
percentage  fluctuation  equals  the  observed  fluctuation  just  mentioned.  The 


adjusted  velocity  wave  components  so  obtained  would  give  us  the  vertical 
and  horizontal  winds. 

For  part  three  we  investigate  the  small  scale  spatial  temporal  airglcw 
fluctuations  produced  by  gravity  waves.  For  this  report  we  shall  only 
consider  a  simple  two-layer  model  plus  a  rigid  ground  boundary .  We  shall 
also  be  exclusively  concerned  with  lew-altitude  sources  which  form  an 
important  contribution  to  the  medium-scale  T.I.D.  's.  The  two-layer  model 
provides  much  of  the  basic  physics  without  any  extraneous  complications.  A 
much  more  detailed  treatment  using  the  OQSPAR  atmospheric  model  plus 
horizontal  wind  profiles  as  well  as  careful  consideration  of  viscosity, 
non-linearity  and  dissipation  is  being  worked  cut.  So  far,  the  two  results 
differ  mostly  in  qualitative  details. 

The  basic  physics  rests  on  the  fact  that  a  gravity  wave  propagating 
upwards  would  reach  height  levels  where  the  frequencies  of  its  higher 
Fourier  components  (within  the  range  of  the  Brunt  frequencies)  become  equal 
to  the  local  Brunt  frequencies  causing  them  to  be  reflected  earthward, 
leaving  only  the  Fourier  cotponents  below  the  Brunt-frequency  range  to 
continue  their  upward  propagation.  The  variation  of  the  Brunt  frequency 
curve  with  height  would  then  serve  as  a  filter  which  filters  out  the  higher 
frequency  components  and  allowing  only  the  longer  period  Fourier  components 
to  propagate  into  the  ionospheric  regions.  In  the  meantime,  the  higher 
frequency  reflected  cotponents  would,  with  the  help  of  the  ground  boundary, 
form  a  complete  set  of  fully  guided  modes  for  each  frequency.  We  shall 


investigate  the  frequency  dispersion  of  the  spectrums  of  fully  guided  modes 
within  the  context  of  a  two-layer  model  as  part  three  of  this  report. 


(II)  Observations  &  Data  Analyses  for  1984  MAPSTAR  Program 

(1)  Phase  Correlation 

In  general  ,  we  assume  that  the  gravity  wave  is  carried  by  the  major 
atmospheric  constituents.  Ihe  response  of  the  minor  constituents  (which 
produce  the  airglow)  through  collisions  with  the  major  constituents  depends 
on  their  vertical  concentration  gradients  relative  to  the  major 
constituents.  If  the  minor  constituent  has  a  layered  structure  such  as  0 
which  produce  the  5577  OI  or  H  and  03  which  produce  the  OH  emissions,  then 
there  is  a  180#  phase  change  between  the  response  at  the  bottom  side  and 
the  top  side  of  the  layer  with  a  continuous  variation  in  the  phase  of  the 
response  in  between  (Thome,  [1968],  Testud  &  Francois  (1971),  KLostermeyer 
(1972  a,b)  Porter  and  Tuan  (1974) ,  Dudis  and  Reber  (1976)  Chiu  and  Ching 
(1978),  Hatfield,  Tuan,  and  Silverman  (1981)).  In  this  section  we  shall 
consider  the  particular  case  of  5577  01  and  7900  A  OH  both  observed  at 
Jelm,  Wyoming.  Since  what  we  observe  is  integrated  along  the  line  of 
sight,  the  actual  phase  of  the  integrated  response  depends  therefore  on  the 
relative  importance  of  the  local  response  at  different  points  along  the 
line  of  sight.  The  situation  is  further  complicated  by  the  fact  that  we 
are  not  in  general  dealing  with  gravity  waves  with  plane  wave  fronts  and 
that  the  source  for  such  waves  may  be  either  very  close  or  some 
considerable  distance  away.  As  a  result,  we  should  not  expect  any  definite 


phase  relationship  between  any  two  different  types  of  uirglcw  for  different 
sources.  We  do  expect  definite  phase  relationship  to  exist  for  the  very 
distant  sources  and  the  long-period  gravity  waves  which  produce  the  large- 


i 
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scale  T.I.D. 's.  A  comparison  of  the  density  variations  in  H  and  O3  with 
the  density  variation  of  the  major  atmospheric  constituents  reveal  the  180# 
phase  change  between  the  bottom  and  top  sides  of  the  H  and  O3  layers 


(Hatfield  et  al  (1981)). 


Fig.  la  shews  the  5577  01  superimposed  over  the  simultaneous  OH 
observations  for  day  176  (June  23,  1984) .  In  this  particular  case  there 
does  appear  to  be  considerable  correlation.  The  variations  appear  to  be  in 
phase.  Since  both  the  ambient  5577  OI  and  the  7900  OH  have  a  layered 
structure,  this  result  would  have  been  expected  for  a  relatively  distant 
source. 


Fig.  lb  shows  the  phase  correlation  by  a  plot  of  the  cross  variance  as 
a  function  of  the  lag  time.  The  positive  maxima  centers  at  s  =  -20  which 
indicates  that  the  5577  01  and  the  7900  OH  are  approximately  out  of  phase 
by  about  20  min.  Both  sums  dropped  off  to  zero  at  about  60  min.  to  either 
side  of  the  center  of  the  maxima  shewing  a  fundamental  4  hour  period  which 
agrees  with  the  crossed  pewer  spectrum  of  the  5577  01  and  the  7900  OH  (dark 


line)  in  Fig.  lc.  The  latter  shews  a  very  high  peak  at  frequency 
f  =  .0038  min.-1  which  also  corresponds  to  approximately  4  hours. 


Fig.  2a  shews  a  plot  of  the  same  two  spectral  lines  for  day  180 
(June  27,  1984) .  Here  we  should  remark  that  even  though  we  have  not  chosen 
this  particular  day  for  detailed  analysis  (there  is  no  radiometric  data  for 
this  day)  we  have  decided  to  show  it  as  the  most  active  day  (at  least  at 
Jelm,  Wyoming)  in  the  entire  campaign.  This  is  evident  both  from  the 
observed  data  and  from  the  computed  percentage  r.m.s.  fluctuations.  There 
does  not  appear  to  be  any  obvious  phase  relationship.  An  additional 
possibility  exists  where  the  relative  strengths  of  different  sources  vary 
throughout  the  night  so  that  phase  correlation  can  also  vary.  The  two 
curves  can  be  in  phase  for  part  of  the  night  and  out  of  phase  for  the 
remainder  of  the  night. 

Fig.  2b  shews  the  plot  of  the  cross  covariance  as  a  function  of  lag 
time  for  June  27  (1984) .  The  negative  minima  which  centers  again  at 
s  =  -15  and  indicates  that  the  5577  01  and  the  7900  OH  are  out  of  phase  by 
about  15  min.  from  purely  180#  out  of  phase.  The  cross  covariance  goes  to 
zero  at  about  55  min.  from  the  center,  giving  us  a  fundamental  period  of 
about  220  min.  The  crossed  power  spectrum  (Fig.  2c)  shews  a  peak  which 
centers  at  about  200-1  min.  which  is  consistent  with  the  crossed-oovariance 
curve. 

Fig.  3a  shews  the  plot  of  the  two  spectral  lines  for  day  182  (June  29, 
1984) .  This  was  chosen  as  the  most  active  day,  although  its  r.m.s. 
fluctuations  is  only  the  second  largest  in  the  campaign.  Once  again,  the 
phase  correlation  appears  sporadic.  A  definite  correlation  appears  to 


exist  beyond  the  3rd  hour  of  the  observation.  Before  that  time,  one  can 
perhaps  make  a  case  for  some  anti-correlation  (i.e.  180#  out  of  phase) . 


|  Since  the  phase  correlation  varies  from  approximately  180#  out  of 

I 

I  phase  to  being  approximately  in  phase,  we  should  use  a  moving  window  of 

!  about  2  hour  duration,  or  divide  the  night  up  into  several  parts.  With  a 

i 

!  single  fixed  window  for  the  entire  night,  the  cross  covariance  (Fig.  3b) 

does  not  shew  much  in  the  way  of  correlation.  In  general,  for  both  of  the 
very  active  days  (June  27  &  29) ,  it  appears  that  phase  correlation  for  the 

! 

I 

i  small  period  fluctuations  appear  to  exist.  For  instance,  on  June  27,  (Fig. 

2a)  considerable  correlation  exists  for  the  15  min.  fluctuations  at  67.5 
min. ,  77  min. ,  and  99  min.  A  phase  shift  begins  to  occur  at  about  2  hours 

I  and  a  sub-peak  in  the  OH  appear  to  be  shifted  by  about  7.5  min.  from  the 

! 

i 

corresponding  peak  for  01.  The  phase  shift  appears  to  be  further  increased 
at  157.5  min.  where  the  peak  for  OH  appears  to  be  shifted  by  about  14  min. 
For  June  29  (Fig.  3a)  the  matching  of  subpeaks  before  3  hours  appear  to  be 
more  difficult,  although  the  OH  and  01  appear  to  be  in  phase  after  3  hours. 

(2)  Frequency  Correlation 

Unlike  phase  correlation  which  depends  on  the  many  factors  mentioned 
in  the  last  section,  the  frequency  correlation  should  be  more  readily 
established,  since  it  is  the  independent  variable  and  the  response  of  the 
minor  atmospheric  constituents  which  produce  the  airglow  must  necessarily 
share  some  common  frequency  spectrum  with  the  atmospheric  disturbance  in 
the  major  constituents. 


As  already  mentioned  in  the  Introduction,  although  we  have  used  four 
different  kinds  of  pcwer-spectrum  analyses,  we  shall  rely  primarily  can  the 
maximum  entropy  method.  The  following  procedure  will  be  adopted  for  each 
set  of  observations:  (1)  We  eliminate  the  unwanted  window  together  with 
its  possible  leakage  produced  by  its  finite  data  base.  Here,  although  a 
number  of  papers  exist  which  deal  with  the  existence  of  the  Window  (Beamish 
(1977)  Gasser  (1974)  Martin  (1979),  Wahbe  (1968),  Beamish  and  Priestly 
(1980),  Ulrich  (1972),  Ulrych  et  al  (1973),  Ulrych  and  Bishop  (1975))  for  a 
maximum  entropy  analysis  (it  certainly  does  occur  for  the  Fast  Fourier 
method) ,  we  have  found  that  it  does  exist  even  if  there  is  relatively 
little  leakage.  To  be  safe,  we  also  discard  the  higher  harmonics 
corresponding  to  each  window  and  shall  refer  to  both  as  "background 
harmonics".  In  the  event  an  exceptional  prominent  peak  occurs  close  enough 
to  one  of  the  background  harmonics  so  as  to  be  indistinguishable,  we  check 
with  the  observed  data  to  see  if  such  a  peak  oscillation  may  indeed  occur. 
If  the  check  is  positive,  we  identify  it  as  a  possible  real  oscillation: 

(2)  We  pick  out  all  the  other  possible  series  whose  fundamental  is  not  a 
member  of  the  background  series.  The  difficulty  lies  in  the  higher 
frequencies  where  the  members  of  the  different  series  begin  to  merge  and  a 
member  of  one  series  may  also  appear  to  be  a  different  member  of  another 
series:  (3)  For  the  very  active  and  highly  disturbed  days,  we  check  if  any 
especially  prominent  high  frequency  peaks  with  periods  in  the  range  of  5-12 
minutes  may  occur,  since  highly  disturbed  days  may  produce  instability  with 
this  type  of  periodicity  (Tuan  et  al  (1979),  Schoeberl  et  al,  (1983)): 


(4)  For  Highly  disturbed  days  we  also  search  for  any  possible  subharmonics 
whose  power  peaks  are  related  in  a  universal  way  (Feigenbaum  (1983)) . 

For  June  23,  we  have  simultaneous  observations  of  5577  01  and  7300  A 
OH  available  for  Wyoming  as  well  as  the  5577  01  and  the  Radiometric  data  at 
Colorado.  Fig.  4a,  4b,  and  4c  shew  the  observed  5577  data  and  the  M.E. 
power  spectrum  analysis  for  this  day.  By  replotting  the  higher  frequencies 
on  a  larger  scale  one  can  shew  that  after  discarding  the  background  series 
with  its  fundamental  at  315  min.,  we  obtain  the  following  possible  sets  of 
series  within  the  limits  of  resolution: 

(1)  T  =  70,  18,  6.4,  5.9,  5.1  min. 

(2)  T  =  30,  7.6,  6.1,  5.1 

(3)  T  =  12.9.  6.4 

(4)  T  =  9.4 

The  relatively  prominent  peaks  (underlined)  which  do  not  seem  to  be 
higher  harmonics  of  other  series  are  70,  30,  12.9  and  9.4  min.  Thus,  the 
analysis  appears  to  shew  that  there  are  long  period  oscillations  of  the 
order  of  70  and  30  minutes  with  short  period  oscillations  of  the  order  of 
12.9  and  9.4  minutes.  Not  shewn  are  a  group  of  peaks  from  7  to  5  minutes. 
We  would  like  to  stress  that  all  these  peaks  are  very  much  higher  than  the 
required  2-4  times  greater  than  the  noise  level.  The  two  short-period 
oscillations  of  12.9  min.  and  9.4  min.  are  more  prominent  than  marry  peaks 
immediately  preceding  them  and  do  not  appear  to  belong  to  members  of  eery 
series. 


For  the  7900  A  OH  observed  at  nearly  the  same  time,  we  have  the 
following  possibly  real  series: 

(1)  T  =  37.5.  18.2.  12.0,  7.2  min. 

(2)  T  =  26,  13.3,  9,  5.4 

(3)  T  =  18.2 

(4)  T  =  7.7 

Fig.  5a,  5b,  and  5c  shew  respectively  7900  A  OH  and  the  M.E.  analysis. 
Again,  we  discard  the  background  series.  The  prominent  peaks  (underlined) 
are  37.5,  for  the  large  periods  and  7.7  for  the  short  periods.  Again  the 
7.7  min.  peak  occurs  as  an  isolated  peak  with  many  small  peaks  of  lower 
frequency  preceding  it.  Both  observations  shew  more  peaks  in  the  70-80  and 
30-40  min.  range.  Both  shew  an  18  min.  harmonic.  Both  shew  small  period 
prominence  in  the  12-14  min.  range  as  well  as  the  7-9  min.  range.  The  real 
question  here  is  how  big  a  role  does  the  interference  of  the  background 
series  play  in  either  case.  The  background  fundamental  for  the  5577  A  has 
a  period  of  315  min. ,  while  the  background  fundamental  for  the  7900  A  has  a 
period  of  300  min.  Furthermore,  the  starting  time  is  not  exactly  the  same. 
We  are  presently  examining  the  data  with  the  same  background  fundamental 
and  the  same  starting  time. 


For  the  5577  A  01  at  Bootlake,  Colorado,  we  obtain  the  following 
possible  series: 

Tx  =  21,5,  10,  6JL  5^6 
T2  =  15.8,  7.9,  5^6 


T3  =  13.1,  6*8,  4*6 

T4  =  12,  5*6 

The  prominent  peaks  (Fig.  6a,  6b,  6c)  appear  quite  strongly  for  the 
shorter  period  oscillations,  particularly  the  6.8  min.  oscillation  which 
stands  out  among  the  rest  of  the  peaks.  A  relatively  broad  peak  occurs  at 
50  min.  (not  listed  above) ,  but  its  height  level  is  merely  comparable  to 
the  6.8  min.  peak. 

For  the  radiometric  observations  in  the  same  day,  the  series  are  given 
by: 

Tx  =  52,  23.9,  16.3,  8*9,  6*2,  4.7 

T2  =  34,  16.3,  8*9,  5.5 

T3  =  7*8 

The  observed  data  and  the  power  spectrum  analysis  are  given  in  Fig. 

7a,  7b,  7c.  The  8.9  and  7.8  min.  peaks  appear  significant  and  are  more 
prominent  than  nearby  peaks.  The  total  useful  data  base  for  both  sets  of 
Colorado  results  is  about  3  hours.  So  we  should  be  concerned  only  with  the 
shorter  period  fluctuations  for  these  data. 

For  the  very  active  day  (Day  182) ,  we  have  simultaneous  observations 
of  5577  A  01  and  7900  A  OH  at  Wyoming,  7300  A  OH  at  New  Mexico  5577  A  01 
and  radicmetry  at  Colorado.  Fig.  8a,  8b  and  8c  shew  5577  A  01  and  the  M.E. 
analysis  for  this  day.  The  length  of  data  base  for  this  day  is  237.5  min. 
Again,  after  discarding  the  background  series  of  frequencies,  we  obtain  the 
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following: 


(1)  T  =  79,  28,  12.8,  11.3,  5^6,  5.3,  4.4 

(2)  T  =  36.5.  12.8,  5.3 

(3)  T  =  17,  8.2,  5^6 

(4)  T  =  5.6 


Hie  relatively  prominent  peaks  occur  at  79  min.  (which  is  very  close 
to  the  second  harmonic  of  the  background  fundamental) ,  28  min.  and  5.6 
minutes.  Although  it  is  within  our  margin  of  resolution,  the  28  minutes 
peak  can  easily  itself  be  a  fundamental  and  does  not  belong  to  the  first 
series  (series  (1)) .  Unlike  the  active  day  the  very  active  day  shews  a 
surprisingly  prominent  peak  relative  to  other  nearby  peaks  at  5.6  min. 


Fig.  9a,  9b,  and  9c  shew  the  respective  results  for  7900  A  OH  taken  at 
Wyoming.  The  length  of  data  base  is  237.5  min.  and  the  possibly  real  sets 
of  series  are: 


(1)  T  =  158.3,  31.7,  17.6.  11.3,  5.3,  4.1 

(2)  T  =  25 

(3)  T  =  17.6.  9.1? 

(4)  T  =  6.0 


The  interesting  feature  here  is  the  big  peak  at  T  =  158.3  min. ;  almost 
exactly  double  the  79  min.  peak  for  the  5577  A  01.  Since  T  =  158.3  min. 
cannot  be  a  harmonic  of  the  background  fundamental  (237.5  min.) ,  it  lends 


greater  credibility  to  the  possibility  that  the  79  min.  peak  is  real.  One* 
again  we  have  a  prominent  short-period  6  min.  peak. 

Fig.  10a,  10b,  and  10c  shew  the  results  for  the  7300  A  OH  observed  at 
Capilla  Peak.  The  length  of  data  base  is  240  min.  The  possibly  real 
series  are: 

(1)  T  =  96,  25.3,  20,  10.9.  9.8,  7J5,  4^2 

(2)  T  =  40,  20,  13.3,  9.8,  6.8,  5.8 

(3)  T  =  15.5,  7^5 

(4)  T  =  10.9.  5.5 

There  are  a  few  split  peaks.  We  have  assumed  that  a  37  min.  split 
peak  is  more  likely  to  be  about  40  min.  Once  again,  there  are  prominent 
peaks  relative  to  their  immediate  neighbors  at  7.5  min.  and  4.2  min. 

Correlations  between  different  data  taken  at  the  same  location  do 
appear  to  exist,  although  they  do  not  seem  to,  from  the  raw  data  without 
the  power  spectrum  analysis.  Based  on  the  very  limited  information,  the 
very  active  day  (Day  182)  does  appear  to  have  more  prominent  short-period 
oscillations  relative  to  their  neighboring  peaks  as  compared  with  the 
average  active  day  (Day  176) .  Feu:  more  observations  will  have  to  be  made 
before  ary  reed  significance  can  be  attached  to  this  result. 

Fig.  11a,  lib  and  11c  shew  the  results  for  the  5577  A  OI  at  Bootlake, 
Colorado.  The  useful  length  of  data  base  here  is  about  5  hours  and  we  see 


that  like  Wyoming,  there  appears  to  be  considerable  activity.  The 
following  series  may  be  obtained  from  the  pcwer  spectrum  analysis: 

Tx  =  210,  26.3.  8.4.  7.5 

T2  =  57,  11.9,  9A,  iL4,  6.3,  5.2 

T3  =  21.8,  10.4.  5.4 

T4  =  18.0,  6.1 

T5  =  15.4.  7.5 

A  rather  prominent  peak  appears  at  26.3  min.  The  higher  peaks  appear 
much  smaller  by  comparison. 

For  the  radiometric  data  (12a,  12b,  12c)  at  Bootlake,  we  obtain  the 
following  possible  series: 

=  221,  36.9.  20.1,  17.0,  8.5,  5*5,  4*5 
T2  =  60,  20.1,  15.5,  10.2.  8.5,  7*_5,  5^9,  5J> 

T3  =  22.9.  11.5,  7*5,  4*5 

The  length  of  data  base  is  about  5.5  hours.  We  have  an  especially 
prominent  peak  at  10.2  min. 

In  general,  the  pcwer  spectrum  analysis  for  day  182  (June  29,  1984) 
appears  to  have  less  prominent  small  period  (5-10  min.)  peaks  than  day  176 
This  appears  surprising  in  view  of  the  fact  that  June  29,  by  most  criteria 


is  a  more  active  day.  Later,  we  will  shew  that  the  r.m.s.  fluctuations 
confirm  this  expectation. 


Tables  I,  II,  III,  TV,  V  and  VI  shew  the  common  peak  (peaks  with 
similar  periods)  for  different  ranges  in  the  oscillation  periods.  For  day 
176,  the  most  common  peaks  appear  in  the  short  period  ranges  (Table  V) . 

The  same  is  true  for  day  182,  although  for  this  day  considerable  common 
peaks  appear  in  11-20  min.  range,  (Table  IV) .  For  both  days,  there  exists 
considerable  less  common  peaks  for  the  long  period  fluctuation,  (Table  I  & 
II) .  One  possible  explanation  is  that  the  natural  "bunching"  together  of 
higher  harmonics  at  the  short  period  end.  However,  this  does  not  explain 
the  profuse  number  of  peaks  in  the  5-6  min.  range.  Although  much  more  data 
is  needed  before  any  definite  statement  can  be  made,  we  believe  that  there 
may  be  a  tendency  for  short  period  oscillations  of  the  order  of  4.5-7  min. 
to  occur  with  greater  frequency  than  other  periodicities.  We  may  add  that 
this  type  of  peaks  occur  without  regard  to  location  (i.e.  irrespective  of 
whether  the  observations  were  taken  in  Colorado,  Wyoming  or  New  Mexico) . 
Some  location  bias  seems  to  occur  for  the  longer  period  fluctuations. 

(3)  Comparison  of  Fluctuation  Magnitude 

Just  as  the  case  for  Phase  correlations  (Section  I) ,  the  magnitude  of 
the  fluctuations  also  depends  on  the  vertical  density  gradient  of  the 
unperturbed  minor  atmospheric  constituent  relative  to  the  major  atmospheric 
constituent  (Porter  &  Tuan  (1974),  Dudis  &  Reber  (1976),  Chiu  &  Ching 
(1978)).  In  general,  if  the  vertical  density  gradients  of  the  minor 


constituents  which  produce  the  airglcw  are  opposite  to  the  correspond ing 
gradients  of  the  major  constituents  which  are  presumed  to  carry  the  gravity 
wave  (e.g.  the  bottom  side  of  a  layered  structure  in  the  minor 
constituents)  the  response  tends  to  be  large.  The  percentage  fluctuations 
on  the  top  side  of  the  layered  structure  on  the  other  hand  would  tend  to  be 
of  the  same  order  as  the  major  constituents  carrying  the  gravity  wave. 

In  computing  the  r.m.s.  fluctuation,  one  way  is  to  put  in  a  polynomial 
approximation  to  represent  the  average  "quiet"  background  and  then  oonpute 
the  average  r.m.s.  deviations  from  this  "background".  Rather  than  making 
this  assumption,  we  have  decided  to  compute  the  percentage  r.m.s. 
fluctuation  about  the  averages  over  60,  30,  15,  10  and  5  min.  intervals. 

The  fluctuations  computed  over  the  5  min.  interval  averages  must  only 
include  the  short  period  oscillation.  In  this  way  we  do  not  have  to 
specifically  assume  some  kind  of  polynomial  "background". 


Table  VII  shows  the  percentage  r.m.s.  fluctuation  for  day  176  and  182. 
The  checks  indicate  the  day  with  clearly  the  greater  fluctuation,  whilst 
the  checks  with  crosses  indicate  the  day  with  the  slightly  to  somewhat 
greater  fluctuation.  A  glance  at  the  Table  would  shew  that  day  182  (June 
29)  is  clearly  the  more  active  day.  A  second  glance  would  shew  that  whilst 
June  29  wee  clearly  the  more  disturbed  day  in  terms  of  long-period 
fluctuation,  June  23  appears  to  have  a  very  slight  edge  for  the  shorter- 
period  fluctuations .  Such  a  conclusion  is  consistent  with  the  power 
spectrum  analyses  which  seemed  to  shew  relatively  l  ore  prominent  short 
period  peaks  for  June  23.  The  last  row  of  Table  VII  shews  the  difference 
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in  percentage  fluctuation  between  the  2  days  expressed  as  a  percentage  of 
the  June  29  fluctuation  percentage.  With  exceptions  in  one  radiometric 
data  at  Colorado  and  one  01  data  at  Wyoming,  the  difference  in  percentage 
fluctuation  is  clearly  less  for  the  shorter  10  to  5  min.  periods.  This 
means  that  the  June  23  data  shewed  relatively  greater  short-period 
fluctuations. 

Table  VIII  shows  the  total  r.m.s.  percentage  fluctuation.  Ihe 
averaging  was  made  over  the  entire  data  base.  The  result  shews  a  very  high 
Colorado  radiometric  data  fluctuation.  Since  this  averaging  includes  the 
"quiet”  background  variations,  one  cannot  make  any  meaningful  more  detailed 
analysis. 

Table  IX  shews  the  percentage  r.m.s.  averages  of  all  different  data 
for  each  day  as  well  as  the  percentage  of  the  shorter  period  oscillations 
to  the  60  min.  oscillation.  Clearly  again  June  23  has  more  significant 
short-period  oscillations. 

In  Table  X,  we  present  a  comparison  of  the  percentage  r.m.s. 
fluctuations  for  six  days  (from  June  21  to  June  29  inclusive)  of  the 
Wyoming  5577  01  and  7900  OH  data.  As  already  mentioned  in  the 
Introduction,  day  180  appeared  to  be  very  active  in  Wyoming  both  from  the 
observed  data  and  from  the  ccnputed  r.m.s.  fluctuations.  Clearly,  the 
column  for  June  27,  Table  X  shews  that  the  percentage  r.m.s.  fluctuation 
for  this  day  is  consistently  the  largest.  The  day  we  chose  (June  29)  as 
the  most  active  is  arguably  the  second  most  active. 


The  most  interesting  consistent  result  is  that  irrespective  of  whether 
the  averaging  was  over  60,  30,  15,  10  or  5  min.  intervals,  the  magnitude  of 
the  percentage  r.m.s.  fluctuation  is  consistently  greater  for  the  5577  01 
as  compared  with  the  7900  OH.  There  is  one  exception  for  the  30  min. 
interval  average  on  June  23  where  the  percentage  r.m.s.  fluctuation  was 
slightly  in  favor  of  the  7900  OH.  Assuming  that  there  are  no  systematic 
errors  or  instrumental  bias,  we  may  conclude  that  for  the  5577  OI,  the 
bottom  side  of  the  atonic  oxygen  profile  has  a  sharper  vertical  gradient  as 
compared  with  the  constituents  which  produce  the  7900  OH,  and  that  on  the 
whole  the  5577  01  profile  is  primarily  concentrated  over  a  narrower  layer 
than  say  the  7900  OH. 

A  similar  statement  may  be  made  about  the  radiometric  data  which  seems 
to  show  from  Table  VII  a  consistent  greater  percentage  fluctuation  than 
either  the  5577  OI  or  the  7900  OH. 

(Ill)  Extraction  of  Velocity  Field  Profiles  From  Airalcw  Observations 
(1)  FFT  Decomposition  of  Observed  Data 

The  observed  optical  emissions  data  may  be  Fourier  decomposed  in  the 
following  way: 


Q(t)  =  Qo(t)  +  AQ(t) 

=  SQ  +  Z  (a^coswfct  +  bj^sinw^t) 


(1) 


where 


Qq  =  unperturbed  background  trend 
AQ  =  fluctuations  produced  by  gravity  waves 
S0  =  time  averaged  unperturbed  trend  (quite  close  to  Qq) 
ak»  ^k  =  Foxier  coefficients  in  the  deconposition  of  AQ(t) 
=  frequency  of  the  kth  harmonic  in  the  deconposition. 

We  may  hence  write, 


AQ  (t , 


a^cosj^t  +  b^sirw^t 


(2) 


where  AQ(t,cjjc)  is  the  observed  fluctuation  for  each  frequency,  Since 

most  available  gravity-wave  models  are  essentially  monochromatic  models,  we 
use  a  model  for  each  frequency,  u^. 


To  determine  the  gravity  wave  amplitude  Which  will  give  the  observed 
fluctuation  given  by  (2) ,  we  consider  the  gravity-wave  induced  fluctuations 
in  airglcw.  To  be  specific,  we  consider  observation  along  the  zenith  for 
the  special  case  of  OH  emissions.  We  consider  the  following  reactions: 


H  +  03  -  OH*  +  02 
OH*  -*  OH  4-  hi/ 


(3) 

(4) 


If  k  is  the  rate  coefficient  for  (3)  and  £  is  the  efficiency  of  photon 
emission  per  OH  molecule,  we  may  write: 
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Where  I  is  the  total  volume  emission  rate.  We  may  write. 


I  =  I0  +  Al 


Where  I0  is  the  unperturbed  volume  emission  rate  and  A I  is  the 


perturbed  emission  rate  produced  by  the  gravity  wave. 


AI  =  fZko[03]o[H]o 


Ak  *jo3_]  MM 

ko  [o3]0  [H]0 


where  a[03],  a[H]  and  Ak  are  the  fluctuations  produced  by  the  gravity 


waves.  In  (7)  we  have  neglected  the  non-linear  terms  for  simplicity. 


Throughout  this  report  we  will  use  not  only  linearized  gravity  wave  models, 


but  also  linearized  response. 


Through  assuming  a  horizontally  stratified  atmosphere, 
a[03]  =  A[03]ei “  ktfO 

A  [H]  =  AfHJei^t  ~  kXx) 

Ak  =  Ake^(wt  “  kxx) 


where  a[03],  a[H]  and  Ak  on  the  right-hand  side  depend  only  on  z. 


Substituting  (8)  into  (7)  and  integrating  over  z  (the  vertical  column) 


we  obtain. 
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If  we  make  the  observation  at  x  =  Xq  and  take  the  real  part  of  (9) ,  we 


obtain, 


S  /  f  Wo^o 


.  4K  +  ±£fs2  +  urn 1  & 

ko  l°3!o  [H1o 


x  (cosw^tcosk xxq  +  siru^tsink^?^) 

s  S  ?  ko[°3>otI,lodz 


Again,  using  linearized  response  and  neglecting  chemistry,  we  may 


write. 


V  . 
phx 


AU  + 


.  f  ,  dN 

i  I  1 - 3  +  9-  . 
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where  [x]  is  an  arbitrary  species.  Here  we  also  note  that  (11)  is  real. 


since  aw  is  pure  imaginary.  The  rate  coefficient  can  be  expressed  in  terms 


of  temperature,  (Jackman  &  McFeters  (1985))  and  is  given  by: 


k  =  Ae 


from  (12)  we  obtain 


■(tils) 
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and  it  is  easy  to  shew  that 


AT  ,  _ .  fa)  A  .  i  “b2  A 
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Thus  all  the  ratios  in  the  bracket  of  (7)  can  be  expressed  in  terms  of 
a  linear  combination  of  au  and  aw. 


(2)  Matching  of  do 


ironochromatic  data  with  monochromatic 


avi tv-wave  models 


If  we  define  a  standard  velocity  field  by, 


AU(Z,  Ufc)  =  QfcAUg  (z,  uk) 
AW  (Z,  wfc)  =  OfcAWg  (z,  Wk) 


We  may  then  adjust  aj^  so  that 


AQ(t,  o,K) 


=  akRe 


where  AQT(t,  new  contains  AUg  and  Aws  only.  Fran  comparing  (16)  with 
(10)  and  (2) ,  we  obtain 


f  S  j 

a*  =  aK4JsooslVio  [  ^  j 
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hfc  -  l  fl  ) 


In  (17)  we  know  a^  and  b^  from  the  Fourier  decomposition  of  observed 
data  and  we  also  know  kx  if  we  knew  the  horizontal  phase  velocity.  Thus, 
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Once  0^  is  determined,  we  can  determine  the  velocity  fields  from 
equation  (15) . 


To  get  a  feel  for  hew  this  can  be  used,  we  shall  assume  that  the 
observations  are  taken  along  the  zenith.  We  will  use  the  starriard 
atmospheric  models  of  Yu  et  al  (1980)  which  is  based  on  OQSPAR.  We  will 
work  with  the  7900  OH  for  our  observed  data  since  for  this  data  the 
variation  in  Qo(t)  for  June  23  and  29  of  1984  is  quite  small  and  may  be 
replaced  by  the  time-averaged  constant,  SQ,  given  by  equation  (1) .  The 
results  for  June  23,  1984  at  periods  of  150  min.  and  37.5  min.  are  listed 
in  Tables  XI  (a)  and  XI  (b) .  Hie  Table  lists  the  velocity  fields  at  90  km 
for  different  horizontal  phase  velocities.  Ibis  is  necessary  since  we  did 
not  knew  the  observed  values  for  the  horizontal  phase  velocities  which  must 
hence  be  treated  as  a  variable  parameter.  These  are  two  particular  and 
relatively  prominent  Fourier  components. 


The  results  for  June  29,  1984  at  periods  of  70  &  19.1  min.  are  given 
in  Table  XI  (c)  and  XI  (d) .  The  velocity  fields  as  computed  from  (16) ,  (9) 
and  (2)  at  90  km  are  listed.  To  construct  the  actual  non-nonochromatic 
velocity  profile,  we  need  to  repeat  the  process  for  every  period  (including 
those  that  do  not  have  a  power  peak) .  The  present  results  appear  to  be 
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rather  low,  but  this  is  expected  since  we  have  only  confuted  two  particular 
Fourier  components  for  each  day. 


In  addition  to  all  the  above  problems,  we  have  used  a  model  (Hatfield 
et  al  (1981))  which  contains  only  the  dynamics  for  our  present  analysis  and 
have  not  used  either  chemistry  or  the  background  winds.  Whilst  recently 
Schubert  &  Walterscheid  (1986)  have  included  more  chemistry,  their  dynamics 
is  based  on  simple  perturbation  with  a  Hines  gravity  model  and  no  wind 
effects  at  all.  So  to  accomplish  the  task  properly,  we  need  to  include  not 
only  all  the  chemistry,  but  also  the  proper  non-perturbative  dynamics  with 
the  inclusion  of  background  winds.  This  is  in  fact  one  of  the  tasks  we 
will  be  working  on  for  our  next  project. 

(IV)  A  Preliminary  Study  on  the  Behavior  of  Short  Period  Gravity  Waves 

(1)  Preliminary  Discussion 

Whilst  a  large  body  of  literature  exists  on  the  guidance  of  gravity 
waves,  the  guidance  mechanisms  for  the  internal  (rather  than  surface) 
gravity  waves  have  been  principally  attributed  to  either  the  structural 
guidance  produced  by  unevenness  in  the  background  atmosphere,  or  the 
unevenness  in  the  background  wind  profile,  (Hines  &  Reddy  [1967],  Reddy, 
[1969]) ,  or  the  "dissipative"  guidance  produced  by  unevenness  in  the 
background  dissipation  (Tuan  [1976],  Richmond  [1978],  Yu,  et  al.  [1980] 

Tuan  and  Tadic  [1982]) .  The  last  guidance  mechanism  has  been  specifically 
investigated  by  Richmond  (1978)  in  a  three-layer  model. 


In  the  present  report  we  wish  to  focus  on  a  particular  guidance 
mechanism  which  can  exist  only  for  short-period  gravity  waves  with  periods 
in  the  range  of  the  Brunt  period  (about  5  to  15  minutes) ,  and  which,  so  far 
as  we  are  aware,  has  not  been  systematically  investigated.  The  basic  idea 
of  the  mechanism  is  that  the  Brunt  period  (the  period  for  natural 
atmospheric  oscillation)  increases  almost  roonotonically  with  height  above 
105  km  (Fig.  13) .  A  gravity  wave  with  frequency  w  from  low  altitude  will 
eventually  reach  an  altitude  zc  =  zQ  (wb)  where  w  =  wb.  Since  for  z  >  zQ, 
w  >  gravity-wave  propagation  becomes  inpossible,  the  wave  energy  will  be 
fully  reflected  earthward  to  produce  together  with  reflections  from  the 
ground  boundary  a  series  of  fully  guided  modes.  Since  most  low  altitude 
gravity  waves  have  horizontal  phase  velocity  Vp^  less  than  the  local  speed 
of  sound,  c  (Battaner  &  Molina  [1980]),  vertical  propagation  in  the 
acoustic  mode  for  z  >  zQ  is  ruled  out  by  the  Hines'  dispersion  relation  for 
Vphx  <  c. 

As  a  preliminary  investigation  to  illustrate  the  feasibility  of  this 
particular  mechanism  we  shall  consider  a  simple  two-layer  model  (Wang  & 
Tuan  (1985))  with  a  rigid  surface  ground  boundary  (Fig.  13).  The 
simultaneous  presence  of  other  guidance  mechanisms  such  as  the  wind 
profile,  the  dissipation,  or  the.  use  of  a  more  realistic  atmosphere  will  be 
postponed  for  a  full  length  paper.  The  present  preliminary  treatment  is 
intended  to  explore  only  the  basic  idea  of  this  ducting  mechanism. 

In  view  of  the  fact  that  the  guidance  mechanism  under  discussion  is 
always  in  effect  whereas  some  of  the  other  guidance  mechanisms,  such  as  the 


wind  profile,  may  not  always  be  in  effect  (the  wind  may,  for  instance,  be 
perpendicular  or  at  a  large  angle  to  the  direction  of  wave  propagation) , 
one  would  expect  that  in  general  simultaneous  observations  of  fluctuations 
belcw  and  above  the  Brunt  period  curve  (Fig.  13)  (e.g.  the  simultaneous 
observations  of  the  OH  or  5577  01  and  the  6300^  red  line  or  the  P.R.D. 
radar,  etc.)  should  reveal  that  short-period  fluctuations  in  the  lower 
region  would  tend  to  occur  more  frequently  than  the  upper  region.  This 
guidance  mechanism  may  then  be  considered  as  a  filter  which  filters  out  the 
lew  altitude  short-period  gravity  wave  components,  leaving  only  the  longer 
period  components  to  penetrate  to  higher  altitudes. 

If  a  more  detailed  calculation  using  a  more  realistic  atmosphere 
should  shew  that  the  suggested  filtering  process  is  significant,  then 
simultaneous  experimental  observations  above  and  below  the  Brunt  frequency 
curve  may  verify  this  process.  We  may  add  that  short  period  oscillations 
can  still  occur  at  hicfli  altitudes  (ewing,  for  instance,  to  high  altitude 
sources  or  to  Doppler-shift  of  longer-period  waves  by  winds,  etc.)  The 
removal  of  these  short-period  oscillations  at  say  the  63O0A  peak  by  high 
altitude  dissipation,  either  through  viscosity  or  heat  conduction,  is 
unimportant,  since  at  that  altitude,  the  dissipation  is  either 
insignificant  (for  Vp^  >  120  m  sec-1),  or  the  wave  amplitude  is  already 
too  small  (because  of  ducting)  to  be  of  interest. 

2.  Method 

We  approximate  the  variation  of  the  Brunt  period  in  Fig.  13  by  two 
layers.  The  atmosphere  is  assumed  to  be  uniform  and  isothermal  in  each 


layer.  The  atmospheric  parameters  for  each  layer  are  taken  to  be  the  same 
as  those  adopted  by  Thane  (1968) ,  although  we  differ  from  Theme  by  adding  a 
rigid  surface  ground  boundary.  The  Brunt  periods  for  these  two  layers  are 
T^i  =  5.4  min  for  the  lcwer  layer  and  =  14.8  min  for  the  top  layer. 

The  advantage  for  such  a  model  is  that  analytic  solutions  are  possible  and 
that  a  Hines’  dispersion  relation  can  be  written  down  for  each  layer.  Thus, 
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where  Vp^  is  the  horizontal  phase  velocity  and  ua  the  acoustic  cut-off 
frequency. 


Quite  obviously,  if  T^i  <  T  <  (such  as  the  eight  minute  period 

wave  in  Fig.  13)  and  Vp^  <  c, 

>  0,  k^  =  (i*)2  <  0  (21) 

For  the  8  min.  wave,  kj  is  a  real  vertical  propagator  whilst  kjj  is 
pure  imaginary.  Within  each  isothermal  layer,  the  hydrodynamic  equations 
can  be  quickly  reduced  to  single  second-order  equations.  The  equations  for 
the  pressure  variation  (normalized  to  constant  amplitude)  are  given  by: 
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where  $  =  wAp/p0  x  a  Ap,  the  actual  pressure  variation, 
z0  =  height  level  which  separates  II  from  I 


The  solutions  can  be  easily  obtained  with  the  following  rigid  surface 
condition  at  z  =  0,  the  finite  energy  requirement  at  z-~°,  plus  the 
interfacial  boundary  condition  at  z  =  zQ; 
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^  -  -J[  j  and  equation  (24)  is  the  rigid  surface 


condition  for  ^  corresponding  to  the  vertical  velocity  AWj  (z=0)  =  0. 
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Applying  (24)  and  (25)  to  solutions  of  (22)  and  (23),  we  obtain 
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Applying  the  interfacial  boundary  conditions  (26)  and  (27)  we  obtain 
the  following  general  dispersion  relation  for  two  uniform  isothermal  layers 
plus  a  rigid  ground  boundary: 
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By  substituting  (19)  and  (20)  into  (30)  we  can  eliminate  kj  and 
kjj  =  it  to  obtain  a  single  transcendental  equation  in  Vp^.  which  can  be 
readily  seen  to  contain  an  infinite  but  discrete  number  of  solutions  for 
each  w,  or  wave  period.  Fig.  14  shows  a  plot  of  some  of  the  lower  discrete 
states  corresponding  to  numerical  solutions  to  the  equation  for  Vpj^  for 
wave  periods  5.5,  6,  7,  8,  9,  10,  11,  12,  13  and  14  minutes.  As  can  be 
seen  from  Fig.  14,  for  sufficiently  low  horizontal  phase  velocities,  the 


discrete  modes  are  so  close  to  each  other  that  they  form  essentially  a 
continuum  for  practical  purposes. 

3.  Remarks  &  Ccmments 

When  a  more  realistic  OOSPAR  atmospheric  model  is  used,  our  ongoing 
investigation  (Wang  and  Tuan  (1985))  would  seem  to  shew  that  whilst  there 
are  important  quantitative  differences,  there  is  little  difference  in  the 
basic  physics  of  the  guidance  mechanisms.  The  results  for  the  first  few 
guided  modes  shew  considerable  resemblance  to  the  two-layer  model  (Fig.  14) 
but  the  locations  of  the  modes  are  different. 

From  Fig.  13,  the  Brunt  period  curve  shews  a  height  level  of  slightly 
below  200  km  at  a  period  of  9  min.  The  CDSPAR-based  calculation  shows  that 
for  horizontal  phase  velocities  ranging  from  96  to  244  m  sec-1,  the  low 
altitude  (<100  km)  waves  with  less  than  9  min  period  become  evanescent  when 
they  reach  200  km  or  less  and  become  negligible  above  220  km.  Using 
airglcw  as  an  exanple,  the  6300A  luminosity  profile  peaks  at  260-290  km  and 
cuts  off  fairly  sharply  below  the  peak  to  about  220  km  and  is  negligible 
below.  Since  the  dissipation  does  not  become  important  until  250  to  300 
km,  where  the  wave  has  long  since  been  filtered  out,  the  suggested  guidance 
mechanism  does  appear  to  be  significant  in  removing  the  short-period 
ccrponents  frcm  the  6300A  red  line  for  lew-altitude  sources.  In  general 
there  are  far  more  lew-altitude  sources  than  high-altitude  (above  100  km) 
sources  and  the  latter  usually  generate  the  large-scale  T.I.D. 's  which  can 


mnrr inrw&nnnruYVw  v 


-V.  VA  S'  V  %"  A"  A  .  A"  A, 


<.‘l 


on  occasion  be  identified  from  their  high  horizontal  phase  velocities 
(-400-700  m  sec"1). 


large  amplitude  waves  can  produce  short-period  oscillations  associated 
with  the  local  Brunt  period  (Tuan  et  al  [1979]).  But  these  oscillations 
vary  with  height  ranging  from  5  min.  at  mesospheric  heights  to  12  min.  at 
ionospheric  altitudes  and  will  not  contribute  to  the  short-period  (<10  min) 
behavior  in  the  upper  region. 


The  presence  of  horizontal  winds  can  alter  the  horizontal  phase 
velocities  of  the  guided  modes  just  as  it  does  for  the  long-period  guided 
modes  (Reddy,  [1969]).  However,  the  probability  that  it  can  have  a 
component  (parallel  to  the  wave  propagation)  large  enough  to  Doppler  shift 
a  9  min.  period  to  less  than  5  min.  at  ionospheric  height  levels  can  be 
estimated  to  be  less  than  35%. 


In  recent  times  there  have  been  considerable  interest  in  short  period 
fluctuations  which  may  be  detected  through  many  different  techniques 
including  mesospheric  airglow  observation  (e.g.  Myrabo  et  al  [1983], 
Petersen  [1979],  Peterson  [1985],  Clairemidi  et  al  [1985]),  or  the  P.R.D. 
radar  system,  (e.g.  Manson  &  Meek  [1980],  or  Manson  al  [1981])  etc).  These 
short-period  fluctuations  are  often  of  the  order  of  minutes  or  lew  tens  of 
minutes,  well  within  the  range  of  the  Brunt-period  curve  (Fig.  13) .  The 
radar  system  can  be  extended  to  accommodate  the  height  range  belcw  and 
above  most  parts  of  the  Brunt-period  curve.  In  the  case  of  airglcw,  since 
the  mesospheric  emissions  are  for  the  most  part  located  below  the 
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Brunt -period  curve,  simultaneous  observations  with  the  6300A  OI  can  readily 
cover  the  range.  The  simultaneous  observations  will  have  to  be  made  at 
less  than  1  minute  intervals.  One  may  then  expect  from  these  discussions 
that  on  an  average  the  frequency  of  occurrence  of  short-period  fluctuations 
in  the  higher  altitude  region  (above  the  Brunt-period  curve)  should  be  less 
than  the  lower  region. 
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TABLE  XI  (d) 

Construction  of  the  horizontal  field  at  90  km 
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